Helioseismology: a fantastic tool to probe the 
interior of the Sun 



Maria Pia Di Mauro 

Teoretisk Astrofysik Center, Bygn. 520, Ny Munkegade, DK 8000 Aarhus C, Denmark 



Abstract. Helioseismology, the study of global solar oscillations, has proved to be an 
extremely powerful tool for the investigation of the internal structure and dynamics of 
the Sun. 

Studies of time changes in frequency observations of solar oscillations from helio- 
seismology experiments on Earth and in space have shown, for example, that the Sun's 
shape varies over solar cycle timescales. 

In particular, far-reaching inferences about the Sun have been obtained by apply- 
ing inversion techniques to observations of frequencies of oscillations. The results, so 
far, have shown that the solar structure is remarkably close to the predictions of the 
standard solar model and, recently, that the near-surface region can be probed with 
sufficiently high spatial resolution as to allow investigations of the equation of state 
and of the solar envelope helium abundance. 

The same helioseismic inversion methods can be applied to the rotational frequency 
splittings to deduce with high accuracy the internal rotation velocity of the Sun, as 
function of radius and latitude. This also allows us to study some global astrophysical 
properties of the Sun, such as the angular momentum, the grativational quadrupole 
moment and the effect of distortion induced on the surface (oblateness). 

The helioseismic approach and what we have learnt from it during the last decades 
about the interior of the Sun are reviewed here. 



1 Introduction 

In the early 60's accurate observations of the photospheric spectrum revealed the 
existence of oscillatory motions, with periods around 5 minutes, on the Sun's sur- 
face [53] , [BS] ■ The observed oscillatory character of the surface was theoretically 
explained by Ulrich [S^] and independently by Leibacher & Stein [12] as due to 
acoustic waves (i.e. p-modes) - generated for some not well known reason in the 
convection zone and maintained by pressure force - trapped in resonant cavities 
between the Sun's surface and an inner turning point, whose depth depends on 
the local speed of sound and frequency. Only few years later more accurate obser- 
vations carried out by Deubner 28 were able to confirm the previous theoretical 
hypothesis about the modal nature of solar oscillations. 

The unprecedented discovery of existence of such phenomenon opened for the 
first time human eyes to the knowledge of the solar interior and formed the basis 
for the development of helioseismology. Like the geoseismology, which studies 
the Earth's interior through the waves produced during the earthquakes, helio- 
seismology study the interior of the Sun through the small oscillations detected 
at the surface. 
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In fact, since each wave, characterized by a specific frequency and wave num- 
ber, propagates through a different region of the Sun, probing the physical prop- 
erties of the crossed medium, like temperature and composition, it is possible to 
deduce the internal stratification and dynamics of the Sun from the spectrum of 
resonant modes. 

Since the first observations, many thousands of modes of oscillation have 
already been identified with great accuracy. The spectrum extends from 0.6 mHz 
to 5.5 mHz and it has a maximum amplitude of about 15 cm/s in velocity. This 
incredible amount of information collected contributed to the success of this 
discipline and has permitted a deep knowledge of the Sun, not imaginable thirty 
years ago. 



2 Theoretical approach to hehoseismology 
2.1 Basic equations of adiabatic oscillations 

A star, like the Sun, is a gaseous sphere in hydrostatic equilibrium, and the 
oscillations are fluid-dynamical phenomena caused by the action of a restoring 
force which arises when the original equilibrium status is perturbed. Such hy- 
drodynamical systems can be described by specifying all the physical quantities 
as functions of the position r and time t. 

To provide a background for the treatment of stellar pulsations, it is useful 
to consider briefly the basic equations of hydrodynamics, which can be derived 
by applying the fundamental principles of conservation of mass, of momentum 
and of energy. We neglect viscosity and magnetic fields and we can also assume 
that gravity is the only acting body force and that the radiation is the dominant 
contribution to the flux of energy. 

The conservation of mass is expressed by the equation of continuity: 

^ = -£> divv , (1) 

where g = g{r, t) is the density and v = dr/dtis the local velocity. The equation 
of motion is: 

dv „ 

g— = -Vp + gV^ , (2) 

where p = p{r,t) is the pressure and * is the gravitational potential which 
satisfies the Poisson's Equation, 

V^* = -AnGg , (3) 

G being the gravitational constant. And finally the energy equation: 

dQ dE d /1\ divF 

dr= d^+^diUJ =^-^' 

where dQ/dt is the rate of heat loss or gain per unit of mass, E is the internal 
energy in unit of mass, s is the rate of energy generation per unit mass, while F 
is the fiux of energy. 
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By using thermodynamic identities the energy equation can be expressed in 
terms of other, more convenient variables, like in |24) : 



dQ 1 f dp Fip dg 



dt g{r3 -I) \dt g dt I ' 
where Fi and are the first and the third adiabatic exponent, defined by: 

where T is the temperature and the derivatives are calculated at constant specific 
entropy. 

The observed amplitude of solar oscillations are very small {6r/RQ ~ 10~^), 
so that the pulsations can be described with accuracy by applying a linear per- 
turbation analysis of the basic equations of hydrodynamic. 

Let us consider a static equilibrium model with pressure po{r)^ density go{r) 
etc. If we consider Eulerian perturbations, the generic physical quantity / can 
be written in the following way: 

fir,t) = Mr) + f'ir,t), (7) 

where fo{r) is the unperturbed term and f'{r,t) is the small perturbation at a 
given spatial point. The small perturbation can also be written in the Lagrangian 
form, by considering a frame following the motion of an element of gas which 
moves from position r to r + 6r: 

Sf{r) = f{r + 5r) - /o(r) = /'(r) + dr ■ V/o . (8) 

Since the typical pulsation period is much smaller than the time required to 
dissipate the thermal energy, we first assume that the adiabatic approximation 
is sufficient to discuss the dynamical characteristics in the interior of the Sun. 
This hyphotesis, which greatly simplifies the treatment of stellar pulsations, is 
not longer verified in the very superficial layers and nonadiabatic effects on the 
frequencies should not be neglected in the study of the surface layers and of the 
pulsation energetics. 

According to these assumptions and by perturbing the four basic equations 
([T|)-(|4]), we obtain the following system of four linear equations in four unknowns 
to study the small oscillations under adiabatic conditions: 

g — — div {ga6r) , (9a) 
dv 

eo-Q^^ -Vp' + e'v* + eoV*' , (9b) 

V^*' ^ -AnGg' , (9c) 

p' + Sr- Vpo - ^ {g + Sr ■ \7go) . (9d) 
Qo 

An extensive and detailed derivation of the equations of adiabatic oscillations is 
provided by Unno et al. in [93j, by Christensen-Dalsgaard & Berthomieu in [T7] 
or by Christensen-Dalsgaard in [TB] . 
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Fig. 1. The spherical polar coordinate system 
2.2 Spherical Harmonic Representation 

The general equations for small oscillations presented above must now be derived 
in the specific case of stars, which are assumed to have a spherically symmetric 
and a time-independent equilibrium structure. 

Let consider a spherical polar coordinates system {r,9,(j)), where r is the 
distance to the centre, 9 is the colatitude, and <j> is the longitude (Fig. [T]). In 
the polar coordinates system a vector field can be written by specifying its 
components in the radial and angular directions. Thus the displacement, for 
example, can be written as: 

dr{r, 9, (j), t) = ^^ctr + S.eae + t (10) 

where a^, o-e and are the unit vectors in the r, 9 and <j) directions respectively. 
By introducing S^h, the horizontal component of the vector, we can also write: 

5r = ^rQr + , (11) 

where S.r = ^rifiO^i'Tt) and = Ch(^:^j0:i) e^i'g the radial and horizontal 
components of the displacement. 

Since the equilibrium state depends only on the radius r, the solutions of the 
linear system (Pa)) - ()9dl) can be obtained in the following form, by separating the 
spatial from the temporal dependence: 

/'(r, 9, 4>, t) = ~f'{r)f{9, cf,) exp(-zwt) , (12) 

where the time dependence has been expressed in terms of an harmonic function, 
characterized by a frequency uj, the amplitude f'{r) is a function of r alone, and 
f{9, (p) describes the angular variation of the solution. 

In the spherical symmetric system, all derivatives with respect to 9 and 
can be expressed in the form of the tangential Laplace operator: 
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Consequently, f{9, (p) can be found as eigenfunction of and it may be chosen 
to be the spherical harmonic Yi"^{d, (p) of degree / and azimuthal order m, which 
indeed satisfies the eigenvalues problem: 

where / and m are integers, such that < m < I, and is the horizontal 
component of the wave number. The spherical harmonics Y™{9,(l)) are defined 
by: 

Vne, </)) = 7V™,iPr(cos0)e^'"'^ , (15) 

where P™ is the associated Legendre polynomial and N^^i is a constant such 
that the following integral over the unit sphere is satisfied: 

27r /'TT 

Yne, 0) sinOdedcP = 5iv5mm' , (16) 

where Sw and 5mm' are Kronecker's deltas, so that the integral is zero \i I ^ I' 
and m ^ m! . 

It follows that the perturbation quantities (jl2p can be written as: 



/'(r, e, ^, t) = y4^/'(r)yr 4>)e-'"' (17) 
and that the displacement vector can be expressed by: 



8r = v47r3fi 



\d9 



smf 



ll"(0,</>)e-*'^* , (18) 



where 3? stands for the real part. 

By substituting the spherical harmonic representations into Eqs. (Pa|) - (l9d[) . 
we obtain the following set of ordinary differential equations, which describes 
the stellar adiabatic oscillations: 



dr \r Apdr/ ' ' gc^ \ijJ^ 

^^,i.^-N^)^. + -^%' + /-f, (19b) 
dr 1 ip dr dr 

1 d f,d^'\ ^^^(p' , e^r,,,\ , ^(^ + 1)^, 



dr V dr 



= -47rG ( + ) + ^ ' ^' , (19c) 

where Si is the Lamb frequency 



^ ^ , (») 



iV is the buoyancy frequency 
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Fig. 2. Illustration of some spherical harmonics. For clarity the polar axis has been 
rotated of 30° into respect the plane of the page. Positive patterns are indicated by 
brighter surfaces while negative patterns are darker 



and c is the speed of the sound in adiabatic conditions, such that under the 
reasonable assumption that the stellar interior can be approximate to an ideal 
gas: 

^2 __ ^ riksT 
Q ~ 



(22) 



where fc^ is the Boltzmann's constant, /i is the mean molecular weight and m„ 
is the atomic mass unit. 

The fourth-order system of ordinary differential equations (I19ap - (|19cl) . to- 
gether with appropriate boundary conditions at centre r = and at the surface 
r = Rq , constitutes an eigenvalue problem, which admits solutions only for par- 
ticular values of the eigenfrequencies uj. The discrete set of solutions, obtained 
for each (Z, m) is labelled with an integer n and describes the so-called spheroidal 
modes. The modes are therefore identified by three quantum numbers, the radial 
order n, which is the number of the nodes of the wave in the radial direction, 
the harmonic degree I, which is the number of nodes on the surface in the di- 
rection of the latitude and the azimuthal order m, which specifies the number 
of the nodes along the longitude on the surface. A few examples of spherical 
harmonics are shown in Fig. [21 It should be noticed, that in the case of spher- 
ical symmetry there are not preferential directions on the sphere, therefore the 
modes show {21 + l)-fold degeneracy in m and both the eigenfrequencies and the 
eigenfunctions do not depend on m. 
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3 Propagation of solar oscillations 



It can be noticed that the set of equations of adiabatic oscillations ()19ap - (|19cp 
can be easily solved, once the boundary conditions are known. However, an 
accurate interpretation of the results requires a more complete description of the 
phenomenon. This can be obtained by an asymptotic analysis of the pulsation 
equations, justified by the fact that the acoustic modes observed in the Sun, 
show fairly high radial order and high degree. 

The asymptotic analysis is usually carried out in the Cowling approximation 
[23] . by neglecting the perturbation ^' of the gravitational potential. In this 
case the oscillations equations p9ap - (|19cl) reduce to a second-order system, that 
according to [2^ it can be written in the following approximate expression: 



1 



If = 



where 

tZ/(r) = c^g^/'^divSr 
The acoustical cut-off frequency loc is defined by: 

2 c2 / dH 
ujj = 7, 1 - 2 



(23) 



(24) 



(25) 



where H — — (d In gi/d?-)^^ is the density scale height. 

In a star in which the main body forces acting are the pressure and the 
gravity, two kind of oscillations can be maintained: the pressure or acoustic 
waves and the internal gravity waves, which form the classes of p modes and 
g modes respectively. The dispersion relation which describes the acoustic and 
gravity waves propagation in a medium, can be derived from Eq. (I23p as: 



(26) 



where fc^ is the radial component of the wave number, which clearly depends on 
the variation of the characteristic frequencies Si, N and Wc with radius. 
The propagation of modes of oscillation requires that > 0: 



oj''-ujl-Sf\l~—]>0. 



It follows that the Eq. (|27p is satisfied in the two domains where: 



2 ^ 2 



and 



(27) 

(28) 
(29) 
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Fig. 3. Propagation diagram of tlie cliaracteristic frequencies A^, u)c and Si, calculated 
for some values of as functions of the fractional radius for a standard solar model. 
The horizontal lines indicate the trapping regions for a g mode with I = 20 and n = 22, 
and two p modes with (Z = 5 , n = 5) and (/ = 20 , n = 7) 



The conditions (j28p and (|29p define the trapping regions of p modes and g modes 
respectively, as illustrated in Fig. |31 Outside these regions the waves are evanes- 
cent and do not show oscillatory character in space and their amplitude decays 
exponentially. 

Detection of g modes would be extremely valuable since they have highest 
amplitudes in the core, and hence their frequencies, if detected, should be very 
sensitive to the structure and rotation of the deeper interior of the Sun. Unfortu- 
nately, although claims for detection of g modes have been made [l^ , we still not 
have any confirmation that they are really excited in the Sun, and the observed 
five-minutes oscillations correspond only to p and f modes. Figure |4] shows a 
set of p modes frequencies obtained in 1996 [75 by the MDI [S? instrument on 
board the SOHO satellite. 

The f modes correspond approximately to surface gravity waves with the 
condition that div(Jr) ~ 0, so that according to Eqs. (Pa)) - (j9dl) . it is possible to 
assume that (5p ~ Jgi ~ 0. The dispersion relation of the f modes is: 

~ 50fc/j , (30) 

where 50 = GM/Rq. Thus, frequencies depend only on the mean density of the 
star, but not on its detailed internal structure. 



Helioseismology: a fantastic tool to probe the interior of the Sun 



9 




1000 







200 



400 



600 



800 



1000 



Fig. 4. A set of p-mode frequencies [75], as function of I, obtained by MDI instrument 
on board of SOHO. Each ridge contains modes with equal values of n 

3.1 Properties of the acoustic modes 

The propagation of p modes in the interior of the Sun can be interpreted very 
simply in geometrical terms, by studying the behaviour of rays of sound, as 
illustrated in Fig. [5l Locally the acoustic modes can be approximated by plane 
sound waves whose dispersion relation is: 



where kr and kh are the radial and horizontal components of the wave vector 
k. This means that the properties of the modes are entirely controlled by the 
variation of the adiabatic sound speed c, which depends on temperature (Eq. l22p . 
From Eq. (|5T|). by using the definition of k^ given in Eq. it follows that: 



At the surface, where c is small, kr is large and hence the wave propagates almost 
vertically. Due to the increase of the sound speed with temperature, kr decreases 
with depth, while kh increases as r decreases, until kr = and the wave travels 
mostly horizontally. This condition is reached at the turning point rt, where: 



uj^ ^c^\kf =c\kl + kl) , 



(31) 





(33) 



At the turning point, the wave is gradually refracted and goes back towards the 
surface. For r < rt, kr is imaginary and the wave decays exponentially. 
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Fig. 5. Propagation of rays of sound in the solar interior in the case of two p modes 
with degrees / = 5 and / — 15. The acoustic waves are reflected at the surface owing to 
the rapid decrease of density, and at the inner turning point rt due to the increase of the 
temperature with depth. Notice that waves with a smaller wavelength corresponding 
to a higher value of the degree I, penetrate less deeply 

It appears clear from the Eq. ([55]) that lower is the harmonic degree I, the 
deeper is located the turning point of the mode (Fig. [5]). Radial acoustic modes 
with I = penetrate to the centre, while the modes of highest harmonic degree 
observed in the Sun {I ~ 1000) are trapped in the outer 0.2% of the solar radius. 

Figure [6] shows eigenfunctions for a selection of p modes with different degree: 
with increasing degree the p modes become confined closer and closer to the 
surface. 

4 Helioseismic investigations 

The frequencies of the solar oscillations depend on the structure of the equilib- 
rium model, predominantly on the local adiabatic speed of sound and in addition 
on the variation of the density and of the adiabatic gradient in the Sun. 

Moreover, the oscillation frequencies have several advantages over all the 
other solar observables: they can be observed with great accuracy and different 
modes probe the characteristics of different layers in the interior of the Sun. 
Thus, accurate observations of the acoustic frequencies, available today from a 
variety of helioseismology experiments on Earth and in space, can be used to 
probe the characteristics and the details of the interior of the Sun. 

The goal of the helioseismology is, in fact, to infer the internal properties of 
the Sun and to understand the physical mechanisms which govern the behaviour 
of our star. This can be pursued by two different complementary strategies. 
The first is the forward approach which consists in comparing the observed 
data with the theoretical frequencies computed for a solar model, following the 
analysis explained in the previous sections. The second is based on the use of 
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Fig. 6. Eigenfunction for p modes with different harmonic degree as function of the 
fractional radius x — v/Rq. Here, the oscillation behaviour is enhanced, by scaling the 
eigenfunctions with the square root of the density and the squared fractional radius 



the observed data to deduce the internal structure and rotation of the Sun by 
means of data inversion. The inverse approach and its results will be extensively 
discussed in the next sections. 

It appears clear that all the helioseismic investigations require the use of a 
solar model resulting from evolution of the structure equations from its formation 
to the present age. The computed models depend on assumptions about the 
physical properties of matter in stars, in particular the equation of state, the 
opacity and the rates of nuclear reactions. 

It is also necessary that the models agree with the known non-seismic prop- 
erties of the Sun: the photospheric radius Rq — (6.9699 ± 0.07) x 10^" cm [T], 
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the observed luminosity Lq — (3.846 ± 0.005) x 10^^ erg [M], the mass 
Mq = (1.989 ± 0.0004) X lO'^'^g as obtained from the study of planetary motion, 
the composition of the photosphere as inferred from meteoritic abundances and 
spectroscopic measurements Z/X = 0.0245 ± 0.0015 [50] , and finally the age 
(4.6 ±0.004) Gy. Furthermore the computation involve some additional hyphote- 
sis and the use of an appropriate theory (e.g. mixing-length) for the treatment 
of the convection, to simplify the theoretical descriptions. 

Here we will show results obtained by using two reference models - Model S - 
by Christensen-Dalsgaard et al. ^18] , which use respectively the OPAL [80] equa- 
tion of state, and the MHD [67] equation of state. The MHD equation of state, 
based on the 'chemical' picture of the plasma, takes into account the effect of ex- 
cited levels of atoms and ions on the properties of plasma and it also considers a 
lowest-order Coulomb coupling term through the Debye-Hiickel approximation. 
The OPAL equation of state, in contrast, is based on a 'physical' description, in 
which nuclei and electrons (free or bound) are the only fundamental constituents 
of the thermodynamic ensemble. 

It is important to point out that much of the uncertainty in a solar model rely 
on the physics which describe the surface, since there are substantial difficulties 
in modelling convective motions and the thermodynamic properties of this region 
as well as in the treatment of non-adiabatic effects on the oscillations. In most 
cases, in fact, the frequencies are calculated in the adiabatic approximation, 
which is certainly inadequate in the near-surface region. 

We will limit here the considerations on the solar modelling, since this is 
not subject of the present review. General presentations and more detailed the- 
ory about computation of standard solar models are described in a number of 
standard texts, e.g. [87]. [2T]. [24 ] . [53 ] . 

4.1 Forward analysis 

A direct way to test a solar model is to consider differences between observed 
frequencies and those calculated for the theoretical model. The aim of this kind 
of investigation is to correct the physics on which is based the solar model in 
such a way to reduce the discrepancies. Among several models, then, we should 
adopt the one which best fits the observed data. 

Historically, one of the main successes of this approach was the spectacular 
overall agreement of the theoretical kh — ui diagram, produced by a standard 
model, with the observed one, showed in the 1988 by Libbrecht [64] . 

Today we have the possibility to handle more accurate helioseismic obser- 
vations. Here, we discuss the results produced by considering helioseismic data 
[75] obtained in 1996 by the MDI instrument on board the SOHO satellite [5^ . 

Figure [7] shows the relative differences between the observed frequencies and 
those calculated by using Model S by Christensen-Dalsgaard et al. [18] as our 
reference model and it employs the OPAL [80] equation of state. The differences 
between the observed and calculated eigenfrequencies at low degree are very 
small and vary slowly with the frequency. However, at high degree the differences 
appear to depend on I and to increase with the frequency. 
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Fig. 7. Relative differences between the observed frequencies [75l obtained by MDI 
instrument on SOHO satellite and the theoretical frequencies computed on the reference 
model, as function of the degree (panel on the left) and of the frequency (panel on the 
right) 

The frequency dependence results in part from uncertainties in the mode 
physics, but also from the real differences between the Sun and its reference 
model. 
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Fig. 8. Relative differences between the observed frequencies p5] and the theoretical 
ones (like in Fig. [T]) scaled by the inertia of the mode and plotted as function of 
frequency 
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The Z-dependence is mainly associated with the variation of the mode inertia, 
since modes with higher I penetrate less deeply and hence have a smaller inertia. 
Thus, high-degree modes are affected more strongly by the near-surface uncer- 
tainties [2] and [49]. The Z-dependence can be isolated in part by considering 
frequency differences scaled by the inertia of the modes or more conveniently 
scaled by Qi that is the inertia of the mode i = (n, I), normalized by the inertia 
of a radial mode of the same frequency (see Fig. [8]). The result is that the major 
inconsistencies, which appear at high frequency, derive from the modelling of 
the surface layers, indicating that the physics applied there, is inadequate for 
describing the relevant phenomena. 

4.2 Solar seismic radius 

The radius of the Sun can be simply determined from measurement of its distance 
and of the angular diameter of the visible disk. But, in practice it is not an easy 
task to distinguish the observed radius from the photospheric radius, defined as 
being the depth where the temperature equals the effective temperature. 

The measurements of the photospheric radius obtained during the past from 
several groups using different instruments have provided results, which appear 
fairly consistent with the standard value, quoted by Allen [1]. 

Only in the 1997, Schou et al. [85] succeed for the first time, in obtaining an 
helioseismic determination of the solar radius by using high-precision measure- 
ments of oscillation frequencies of the f modes of the Sun, obtained from the MDI 
experiment on board the SOHO spacecraft. They determined that the seismic 
radius is about 300 km smaller than the model radius. A similar conclusion was 
reached by Antia [3^ on the basis of analysis of data from the GONG network. 

The helioseismic investigation of the solar radius is based on the principle that 
the frequencies of the f modes of intermediate angular degree depend primarily 
on the gravity and on the variation of density in the region below the surface, 
where the modes propagate. From the asymptotic dispersion relation pO[) . one 

3/2 

can easily deduce that uj cx Rq . Therefore, by applying a variational principle, 
we can obtain a relation between f-mode frequencies, cjj.o — '^'^^i.Oy and the 
correction AR that has to be imposed to the photospheric radius Rq assumed 
for the standard solar model: 



where <> denotes the average weighted by the inverse square of the measure- 
ment errors. 

Dziembowski et al. [3 7) , by analyzing a long time-series of MDI f-mode fre- 
quencies, inferred temporal variation of the solar radius, with the aim of deter- 
mining a possible solar cycle dependence. Their first results (Fig.[9l), covering the 
period from May 1996 to April 1997, show that the maximal relative variation of 
the solar radius during the observed period was about AR/Rq = 6 X 10-6, which 
corresponds to approximately AR = 4 km. Recently, Dziembowski et al. |38j an- 
alyzing a larger set of data spanning a period from mid 1996 to mid 1999, have 




(34) 
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Fig. 9. Relative differences in radius of the Sun inferred from the f-mode frequencies 
and the standard solar model as in |37) 

found that the systematic trend of AR/Rq is not correlated with the magnetic 
activity. 

However, Brown & Christensen-Dalsgaard [TT], by combining photoelectric 
measurements with models of the solar limb- darkening function, determined a 
photospheric radius of 695.508 ± 0.026 Mm. This value appears even smaller 
than the helioseismic one and the reason for this discrepancy is still unknown. 

It is clear that the problem that remains to be clarified is the connection 
between the seismic radius determined from helioseismic measurements and the 
definitions of solar radius as obtained from the other methods. 

5 Helioseismic inversion 

The inverse problem, always associated with the forward approach, involves esti- 
mating some functions to describe the physical properties of the Sun, by solving 
integral equations which appear expressed in terms of the experimental data. In- 
version techniques are well known and applied with success in several branches of 
the physics from geophysics to the radiation theory, as reported in |7D] , [2^ and 
[50] . Applications to the helioseismic data have been studied extensively during 
the last decade and inversion methods and techniques have been reviewed and 
compared by several authors, e.g. [48], [27], [55], [39], [4] and [47] . 

The observed data are related to the physics of the solar structure, in a very 
complicated way, and the main difficulty arises from the fact that the helioseismic 
inversion is an ill-posed problem: 

• For each set of data, there exists an infinite number of solutions 

• The solution is not unique 

• The solution does not depend continuously on the data. 
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In fact, the observed frequencies constitute a finite set of data and the errors 
in the observations prevent the solution from being determined with certainty. 
Thus, an appropriate choice of a suitable technique of inversion is the first im- 
portant strategy to adopt during a helioseismic inverse analysis. 

The first attempts at inversion used analytical methods to solve integral 
equations obtained in first approximation by applying the asymptotic dispersion 
relation of solar frequencies, the so-called Duvall law [34j. This inversion method 
is considered however not very accurate, since the Duvall law represents a rough 
approximation for the low degree modes which are the more appropriate to study 
the Sun's core. For this reason, here we will consider only numerical techniques 
of inversion. 

5.1 Inversion Techniques 

Since most of the fundamental aspects of inversion do not depend on the dimen- 
sions of the space in which the problem is posed, for simplicity here we consider 
the general case of the linear one-dimensional inversion problem, in which the 
measured data di are functionals of a single function, /(r), of the distance to 
the centre r: 



where Rq is the radius of the Sun. The properties of the inversion depend both 
on the mode selection i = (n, I) and on the observational errors e^, which char- 
acterize the mode set (i = 1, . . . , M) to be inverted. 

The observational errors Si in the data, are assumed to be independent and 
Gaussian-distributed with zero mean and variances af. Given a set of data and 
errors, the problem is to determine f{r) by solving the Eq. ([55]) . where IC'^{r), 
the kernels of the integral, are known functions which depend on the quantities 
of the reference model and its eigenf unctions. 

There are two important classes of methods of obtaining estimates of /(r): 
optimally localized averaging method based on the original idea of Backus & 
Gilbert [S], [5] and regularized least-squares fitting method due to Phillips [7T] 
and Tikhonov [91]. Both methods give linear estimates of the function /(r) 
and give results in general agreement, as it was demonstrated by Christensen- 
Dalsgaard et al. [20] and by Sekii [88] . 

Optimally localized averaging (OLA) The localized averaging kernel method 
allows us to solve Eq. (j35|) by estimating a localized weighted average of the un- 
known generic quantity f{r) at selected target radii rg's by means of a linear 
combination of all the data d^: 




(35) 




(36) 



Helioseismology: a fantastic tool to probe the interior of the Sun 



17 



where ai(ro) are the inversion coefficients to be found and 



M 



K{r^,r)=Y,a,{r^)lC'{r) (37) 

i=l 

are the so called 'averaging kernels'. 

Because of the ill-conditioned nature of the inversion problem, it is necessary 
to introduce a regularization procedure. By varying a trade-off parameter 6, we 
look for the coefficients ai(ro) which minimize the propagation of the errors and 
the spread of the kernel: 

/ J(ro,r)if(ro,r)2dr + tan0Va2a2(ro) , (38) 

assuming that 

K{ra,r)<lr = l. (39) 

J{rQ,r) is a weight function that is small near tq and large elsewhere, assumed 
to be: 

J{ro,r)^{r~raf . (40) 
From Eq. we obtain the expression for the inversion coefficients: 

«.(^o) = ^[^», (r-o)+tan0£;,,]-i / JC^ {r)dr , (41) 







where 



2 r 



S^Aro)^ ir~roy[IC\r)IC^{r)]dr , (42) 
Jo 

and the diagonal covariance matrix of the errors has elements: 

p _/o-?a|(''o) fori=j , , 

~ \ for « ^ j ^^"^^ 

By lowering the trade-off parameter it is possible to obtain more localized 
averaging kernels closer to the nominal point r — rg, but this decreases the 
accuracy with which the solution is determined, since the importance of the 
errors increases. Thus, we should choose, among all the possible solutions, an 
optimal compromise between localization and accuracy of the solution. 

The Eq. (|4ip is equivalent to solve the following set of linear equations: 

A(ro)a(ro) = b , (44) 

where Q:(ro) represents the vector for each target radius, whose M elements are 
the coefficients ai(ro); b is the vector which contains the Lagrangian multipliers; 
A(ro) is the M x Al symmetric matrix, whose elements for each tq are = 
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Sijira) + tanOEij. Therefore, the OLA method is very much demanding on 
computational resource, since it requires the inversion of N matrixes of order M 
to determine the solution at N radial points. 

The errors of the solutions are the standard deviations calculated in the 
following way: 



Sfiro) 



■ M 

^a^(ro)crf 

.i=l 



1/2 



(45) 



while the radial spatial resolution is assumed to be the half-width at half- 
maximum of the resolving kernels. 

The same method can be applied in the variant form described by Pijpers 
and Thompson in [73| and [74], known as SOLA method (Subtractive Optimally 
Localized Averaging), making attempts to fit the averaging kernel to a target 
function, usually a Gaussian function G(ro, r), of appropriate width and centered 
at the target radiulos: 

K{ro,r)c^G{ro,r)c^5{ro,r) . (46) 

In this case, the trade-off parameter is rescaled at each target location to keep 
constant the width of the averaging kernels and to obtain more localized resolving 
kernels closer to the nominal concentration point. Therefore, the coefficients are 
determined by minimizing the following: 

(47) 



so that 



where 



•^0 1=1 z=l 



a^{ro) - y^iU,, + tanOE,,)-^ / /C'(r)G(ro, r)dr , (48) 
.Re 

U^J = / [/C*(r)/C^(r)]dr . (49) 







This inversion problem appears equivalent to solve the following set of linear 
equations: 

Wa(ro) = g(ro) , (50) 

where W is a matrix whose elements are wij = Uij -\- tan OEij , and hence does 
not depend on tq; S'(?'o) is the cross-correlation vector of the kernels with the 
target function G(ro,r). Thus, the solutions are obtained by inverting the ma- 
trix W only one time, such that the computational efforts is therefore reduced 
substantially. 



Helioseismology: a fantastic tool to probe the interior of the Sun 



19 



Regularized Least-Squares Fitting (RLS) This method ahows to find a 
solution that is expressed as a linear combination of a chosen set of base functions 
with j = (l,...,iV): 

N 

where fj are constants to be determined. 

We can choose the base functions, for example, being piecewise constant 
functions on a dissection Vj of the interval [0, Rq]: 

'^^(^) = |o elsewhere ^^2) 

so that /(r) = fj on the interval [rj-i, rj]. Other common alternatives are to 
choose (r) as a continuous set of piecewise linear functions or as a set of splines. 

The parameters fj are determined by a least-squares fit to the data. However, 
this procedure needs a regularization procedure to obtain a smooth solution. So, 
basically we can determine the constants by minimizing: 

+ M r^-^/»]'dr, (53) 
Jo 

where /i is a trade-off parameter between resolution and error and is a differ- 
ential operator so that J^f{r) is a suitable weight function that determines the 
relative importance of smoothing in different regions. 

The minimization of the (f53| leads to a set of linear equations which permits 
to determine fj and hence the solutions. 



M 
i=l 



IC'{r)f{r)dr 



6 Inversions for the solar structure 
6.1 The variational principle 

The numerical inversion of data to determine the solar structure is based on the 
use of the variational principle of Chandrasekhar [12]. Thus, the eigenfrequencies 
can be determined by solving an eigenvalue problem, whose expression can be 
obtained directly from the basic equations governing linear adiabatic oscillations 
Eqs. (iniD-(l9d|: 

uj^Sr = J^{Sr) , (54) 

where are the eigenvalues, is a linear operator on the eigenfunctions Sr. 

Although the frequencies of solar oscillations can be known from observations, 
the eigenfunctions cannot be determined experimentally, so Eq. define a 
nonlinear integral equation. 

However, Eq. (|54l) can be linearized around a known reference model, under 
the assumption of hydrostatic equilibrium. This procedure, whose details can be 
found, e.g., in [93', provides a linear integral equation that can be used in an 
inverse procedure to determine the corrections which have to be imposed to the 
reference model in order to obtain the observed oscillation frequencies uji — 2T:ui. 
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6.2 The surface term 

Non-adiabatic effects and other errors in modelling the surface layers, that can 
give rise to frequency shifts, as it was explained in Section 4.1, have to be taken 
into account by including an arbitrary function of frequency Fsurf(j^) in the 
variational formulation, as suggested by Dziembowski et al. in |40) . 

The function -Fsurf('^) must be determined as part of the analysis of the 
frequency differences. It should resemble, in practise, the differences Qi5vi/vi 
plotted in Fig. ([8]). As Fsmiit^) is assumed to be a slowly varying function of 
frequency, it can be expressed as expansion of Legendre polynomials P\{v), 
usually of low degree A. 

In the inversion procedures it is common use to suppress the surface term 
[27] . This is done by constraining the inversion coefficients to satisfy: 

M 

^a^PA(^^^)0^' =0 A = 0,l...yl. (55) 

1=1 

The maximum value of the polynomial degree, yl, used in the expansion is a 
free parameter of the inversion procedure. In practice, we should consider an 
appropriate value of A for any given data set. 



6.3 Inversion for sound-speed and density 

It follows, from the preceding discussion, that the differences in, for example, 
sound speed c and density g between the structure of the Sun and the reference 
model j(? ^ ^qI Q) can be expressed by the following integral equation [40) : 

^ - r A-;.(0^(.)d. + t KHr)'Ar}<ir + + , (56) 

l^i JO C Jq g (^i 

where K^2 g and K^^ ^2 are the kernels. The term Qi has been already introduced 
in Section 4.1. 

Equation (|56p forms the basis for the linearized structure inversion. Unlike 
the case considered in Section 5.1, this linearized inverse problem involves three 
unknown functions: bf? jc?^ 5qI g and Fsuii{v)- However, the number of the un- 
known functions can be reduced to one by adapting the method of the optimally 
localized averages. 

The principle of the inversion, by generalizing the SOLA technique (see Eq. 
I47p . is to form a linear combinations of bvijvi with coefficients oiiiro) chosen to 
minimize: 

/ [/C(ro,r)-5(ro,r)]2dr + /3 / C^iro,r) f{r) dr + J2(^Uro)<yL (57) 
Jo Jo ,_i 



where 



M 

/C(ro,r)=^a,(ro)i^:.^^(r) 

i=l 



(58) 
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are the averaging kernels, while 

M 

C{r^,r) =Y,a,{ro)KlAr) (59) 

i=l 

are the cross-term kernels. The parameter /3 control the balance between the 
contribution from 5q/ g on 8(? jt?; n is the trade-off parameter, determining the 
balance between the demands of well-localized kernels and a small error in the 
solution; /(r) is a suitably increasing function of radius aimed at suppressing 
the surface structure in the cross-term kernel, e.g. we can use /(r) = (1 -|-7-/i?)^. 

Thus, if our goal is to infer the speed of the sound, the coefBcients ai(ro) 
should be chosen such to suppress the contribution from the cross term, to 
localize the averaging kernel near r = tq, to suppress the surface term assuming 
the (j55|) . while limiting the error in the solution, by the use of the two parameters 
/3 and /i. 

Inversion results The first significant results concerning the application of the 
inversion technique to the Sun were obtained in 1985 by Christensen-Dalsgaard 
et al. (19j, who produced the sound speed profile in the interior of the Sun and 
who first determined the location of the base of the convection zone. Since then, 
several efforts have been done for inverting data in order to test the correctness of 
the standard models in view of the improvements accomplished in the description 
of the relevant physics. A significant progress, in particular, has been achieved 
with the inclusion of diffusion of helium and heavy elements at the base of the 
convective zone p]. 
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Fig. 10. The relative squared sound-speed difference between the Sun and the standard 
solar model [IB] as obtained by inversion of MDI/SOHO data [HS] 
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Fig. 11. The relative squared density difference between the Sun and standard solar 
model [18] as obtained by inversion of MDI/SOHO data [83] 



The resulting profiles for the speed of the sound and for the density, which 
are shown here, have been obtained by inversions of high quality helioseismic 
data obtained during 1998 by Schou [83], from SOI-MDI [82] instrument on 
SOHO satellite. This set includes only modes with harmonic degree I < 100. 

The Model S [TH] which employs the OPAL [50] equation of state is used here 
as reference model. 

Figures [TU] and [TT] show the behaviour of the relative squared sound-speed 
and density differences between the Sun and the standard solar model as func- 
tion of the fractional radius. The vertical error bars correspond to the standard 
deviations based on the errors in the mode sets, calculated by Eq. (|45|) . whereas 
the horizontal bars give a measure of the localization of the solution. The results 
indicate the substantial correctness of the standard solar models. In fact, it is 
clear that deviations are extremely small, except below the base of the convec- 
tion zone (0.71 Rq) where the theory fails to correctly describe the turbulent 
convection. 

The structure of the core, however, is still quite uncertain since the few 
modes with lowest harmonic degree that are able to penetrate towards the centre, 
sample the core for a relative short time because of the large sound speeds there. 



6.4 Inversion for equation of state and the solar helium abundance 

The equation of state can be investigated through the first adiabatic exponent Fi , 
the partial logarithmic derivative of pressure with respect to density at constant 
specific entropy, already defined in Eq. ([5]). 

The solar plasma is almost an ideal gas, and the first adiabatic exponent is 
therefore close to 5/3 in most of the interior. It deviates from this value in the 
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zones of hydrogen and helium ionization, near the surface. Therefore, inversions 
of helioseismic data can be used, in particular, to study the equation of state 
and to probe the helium abundance in the solar envelope, as it was proved e.g. 
in [46], [56], [42], [S]. 

An integral equation analogous to Eq. (|56p can be derived to determine 
the behaviour of ((5A/A)int the relative intrinsic difference in Fi, at constant 
pressure p, density g and composition, between the equation of state of the Sun 
and the one of the reference model, as in [8]. 

The kernels for (c^ , g) which appear in Eq. ([56| can be converted to kernels 
for the set (Ji,u,y), where u = p/g and Y the helium abundance. After the 
conversion, Eq. (j56l) can be written as 



where (^-r'i/-ri)int is the difference in Fi that results from the differences in the 
equation of state alone, but not from the resulting change in solar structure. The 
term SY denotes the difference of the helium abundance in the convective zone 
between the Sun and the model. 

According to the the OLA technique of inversion (Section 5.1), the coefficients 
are found by minimizing: 



The parameters /3i and /32 control the contributions of Su/u and 6Y, respec- 
tively, and ^ is a trade-off parameter which controls the effect of data noise. As 
in Eq. (|38)) J(ro,r) is a weight function; /(r) is included to suppress surface 
structure in the first and second cross-term kernels, like in Eq. (j57p . 

Figure [12] from Di Mauro & Christensen-Dalsgaard [30] , shows the resulting 
intrinsic differences in Fi between the Sun and the two available equations of 
state (OPAL and MHD), as obtained by inversion of the data set by Schou [S3], 
which includes only modes with low and intermediate harmonic degree {I < 100). 

As already shown by Basu & Christensen-Dalsgaard [8j, by using only low 
and intermediate-degree modes it is difficult to judge the significance of the 
differences between the two equations of state. Nevertheless, Fig. [T^] confirms 
previous findings by Elliott & Kosovichev [33] that Fi deviate from 5/3 in the 
central core, probably due to relativistic effects. 

The results shown in Fig.jTS] as found in Di Mauro & Christensen-Dalsgaard 
[31] ■ have been carried out by inversion of preliminary helioseismic data by 




(60) 




(61) 
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Fig. 12. The intrinsic difference in the adiabatic exponent Fi between the Sun and the 
OPAL [80] equation of state (filled circle) and the Sun and the MIfD equation of state 
[18) (open triangles) obtained by inversion of a set of data by Schou [S3] , which does 
not include high-degree modes 



Rhodes et al. [7S|, which include high-degree modes {I < 1000), obtained in 
1996 by the MDI instrument on board the SOHO satellite. The set is made 
up of a very large number of data (7480 modes), which makes the computations 
slow and very demanding in terms of computer memory. The precise high-degree 
modes are able to determine variations very near the solar surface, through the 
He II ionization zone and also part of the He I ionization zone, while by using 
only low and intermediate-degree modes (Fig. ll2[) . we cannot determine solutions 
above r ~ 0.967?©. 

From Fig. [T3]we can affirm that, as noticed by Basu et al. [TO], the OPAL 
equation of state is able to describe better the plasma conditions in the interior 
of the Sun below 0.97 Rq. In the upper layers above O.97i?0, the results indicate 
a large discrepancy between the models and the observed Sun, even considering 
higher-order asymptotic terms in Fgnrf ■ Here, the use of very high degree modes 
reveals that the differences between the two equations of state are very small. 
This is in contrast to earlier results by Basu et al. [TO], which found evidence, 
by inverting a set of data with no highest degree modes, that MHD models give 
a more accurate description of the very upper layers than the OPAL models. 

The helium abundance in the solar envelope It is well known that spec- 
troscopic measurements of the photospheric abundance of helium {Yq) in the 
Sun are very uncertain and, before the advent of helioseismology, the only accu- 
rate method to quantify Yq was based on a calibration of solar models, in which 
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Fig. 13. The intrinsic difference in the adiabatic exponent Fi between the Sun and 
the OPAL ^SOj equation of state (filled circle) and the Sun and the MHD ^67j equation 
of state (open triangles) obtained by inversion of a set of data by Rhodes et al. [75] , 
which includes high-degree modes 



the helium abundance has to be adjusted to match the observed solar lumi- 
nosity. The value of helium abundance calibrated with this method is typically 
about 0.27. 

In the 1984, Gough [JB] noted that the strong sensitivity of acoustic modes 
to the variation of the adiabatic exponent in the Hell ionization zone could 
allow also a seismic determination of the helium abundance in the outer layers 
of the Sun. Thus, equations like (|6ip . may be inverted to determine SY, the 
difference between the helium abundance of the Sun and that of the solar model 
in the helium ionization zones [1(7. Since the convection zone is fully mixed, this 
provides a measure of the value of the helium abundance in the solar envelope. 
It is also important to point out that the determination of the solar helium 
abundance, inferred from the inversion of data, is sensitive to the equation of 
state employed in the reference model. 

The first seismic measures of 1© obtained by Christensen-Dalsgaard et al. 
[15] reported values between 0.24 and 0.25, that were significantly less than the 
abundance estimated by the calibration on the standard solar model. Dziem- 
bowski et al. 41 pointed out that the difference was in rough agreement with 
that expected by the effect of gravitational settling of helium and heavy ele- 
ments, as calculated by Cox et al. [25]. So, today settling is contained in all the 
most accurate standard solar models. 

Recently, Di Mauro & Christensen-Dalsgaard [STj have used Eq. (|6T]) to de- 
termine 5Y, by inverting a set of data with high degree acoustic frequencies [75] . 
By using the MHD equation of state, they obtained a value of 0.2426 ± 0.0005, 
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consistent with the earher resuhs by Kosovichev [57] and Richard et ah [75], 
which employed a similar variational technique. By considering the OPAL equa- 
tion of state they obtained a value of 0.2648 ± 0.0004, which is strikingly higher 
than previous values quoted in the literature (~ 0.242 — 0.25) by Basu & Antia 
[7], Kosovichev [57], Richard et al. [77] and Basu et al. [TO] . 

The rather high value obtained for the helium abundance based on the OPAL 
model may be due to the use of a set with very high degree modes. The determi- 
nation of observational frequencies for high-degree modes still suffers from sub- 
stantial difficulties, related to the merging of power into ridges and the proper 
treatment of the leakage matrix (e.g. [83]) and this could cause systematic errors 
in the frequencies. 

7 Dynamics of the Sun 

7.1 Fine structure in the acoustic spectrum of oscillations 

So far, we have considered only oscillations of a spherically symmetric structure, 
but it is well known and easily observed at the photosphere that the Sun is a 
slowly rotating star. 

The rotation breaks the spherical symmetry of the solar structure and splits 
the frequency of each oscillation mode of harmonic degree I into 2/ + 1 compo- 
nents. Multiplets with a fixed n and / are said to exhibit a frequency "splitting" 
defined by: 

= w(n, /, m) — w(n, /, 0) , (62) 

somewhat analogous to the Zeeman effect on the degenerate energy levels of an 
atom. 

The determination of the splittings is often very difficult, so generally, the 
observations have not been applied in terms of individual mode frequencies, 
but rather it is customary to represent the frequency splittings by polynomial 
expansion in terms of the so called a-coefficients, as explained in [84) : 

w(n,/,m) =w(n,Z,0)-H27r ^ aj{n, l)vf\m), (63) 

where T-'j'-' (m) are orthogonal polynomials that can be chosen, for example, 
like Ritzwoller and Lavely in [78j . Because of the symmetry properties of the 
splittings, the solar rotation is described only by the odd coefficients a^, while 
the even coefficients are a measure of the Sun's asphericity. 

7.2 Inversion for solar rotation 

To study the dynamics of the Sun we need to reconsider the derivation of the 
basic oscillation equations (|19ap - (|19cp by including the effect of a velocity field. 
We assume that the rotation is sufficiently slow that the centrifugal force and 
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other effects of second and higher order can be neglected. This treatment allows 
to define a new expression which relates eigenfrequencies with eigenfunction and 
physical quantities, like the Eq. ([54|) . By applying standard perturbation theory 
to the eigenfrequencies, it can be shown that the rotational splittings are related 
to the rotation rate f^{r, 9) inside the Sun by: 

Au;n,i,m= I ° I ^"^''™(r)/2(r,0)rdrd^ (64) 

where 9 is the colatitude and /C"'''™(r) are the mode kernel functions. The 
dependence of the splittings on angular velocity can be used in a 2-dimensional 
inverse problem to probe the dynamics of the Sun. 

The 2-dimensional inverse problem can be simplified by considering that the 
expansion of the splittings in polynomials given in the Eq. ()63p corresponds to 
an expansion of f2{r,9) such that: 

nir,^^)Jf:n,,Mr)^^^ (65) 

where P2j+i (/i) are the Legendre polynomials with fi — cos 9. So, the a-coefficients 
are related to the expansion functions f22j+iif) by: 

27ra2j+i(n, = J /C"' (r)r22j+i (?')dr (66) 

in which the kernels are calculated according to the expressions given in [33] . 

Equation (|66|) constitutes the basis for the 1.5-dimensional inversion. Now, 
the original inverse problem Eq. (|64|) has been decomposed into a series of 1- 
dimensional independent inversions for each a-coefficient to determine the ex- 
pansion functions f22j+i{r), whose combination according to Eq. (j65p leads to: 

r2(r,M) = n,{r) + 4(0^^ + 4(0^^ + (67) 

dfj, dfi 



7.3 Inversion results 

The variation of the Sun's angular velocity with latitude and radius shown here, 
has been determined by Di Mauro et al. 33 by means of a 1.5 dimension SOLA 
helioseismic inversion of more than 30,000 p-mode splitting coefficients. These 
data were obtained from the first set of uninterrupted Doppler images from SOI- 
MDI (on board the SOHO satellite) in 1996 [86 , which yield splittings of great 
accuracy, never obtained in previous sets of data. 

The inferred rotation rate is shown in Fig. [U where the points indicate the 
angular velocity at various depths calculated at the equator, and at latitudes of 
30°, 60° and 75°. 
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In Fig. [15] contours and red-scale indicate isorotation surfaces in a cut of the 
interior of the Sun. The results confirm the previous findings that the latitudi- 
nal differential rotation observed at the surface persists throughout the convec- 
tion zone, while the radiative interior rotates almost rigidly at a rate of about 
A30nHz. At low latitudes the angular velocity, through the largest part of the 
convection zone, decreases with the radius while at high latitudes increases in- 
wards. The near-surface behaviour agrees with the observed surface rotation 
rate. 

The tachocline, the transition layer from latitudinally-dependent rotation 
to nearly independent rotation [S^, is of very considerable dynamical interest. 
Furthermore, it is thought that the global dynamo behaviour, responsible for the 
solar 11 years magnetic cycle, rises from strong toroidal magnetic fields generated 
by rotational shear in this thin region. 

The tachocline appears mostly located in the radiative zone at a pretty sharp 
midpoint near about r — 0.693 i?© according to Corbard et al. [22], and near 
r = 0.695 i?o for Charbonneau et al. [15]. It is also a fairly thin layer, not more 
than 0.05 Rq at the equator. The layer seems to be wider at high latitudes, but 
certainly less than 0.1 Rq [32]. Charbonneau et al. [15,^, have recently confirmed 
that the width of the tachocline appears to change with the latitude, with a 
minimum value at the equator of (0.0039 ± 0.0013) Rq. 

Another interesting dynamical feature occurs near the poles, where unfortu- 
nately it is very difficult to localize the inversion solutions. Fig. [T3] shows the 




r/K. 



Fig. 14. Differential rotation at four latitudes as obtained by a 1.5 dimensional SOLA 
inversion of the SOI-MDI data. The approximate base of the convection zone is indi- 
cated by the dashed line [33] 
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Fig. 15. Rotation rate in the Sun obtained by inversion of MDI data. Colours and 
contours indicate the isorotation surfaces. The white area indicates the region in the 
Sun where the data have no reliable determinations jSSj 

presence at latitude of 75° of a fairly localized region rotating faster than the 
surroundings [86 . It is still not clear, if this feature is somewhat related to the 
applied inversion technique. 

Very recently, Howe et al. [51] have found evidence that the rotation rate near 
the base of the convective envelope shows variations with time, with a period of 
the order of I.Syr at low latitude. Such variations occur above and belove the 
tachocline and appear more pronounced near the equator and at high latitudes. 

To infer accurately the rotation in the deepest interior, it is necessary to in- 
vert a set of data which includes accurate splittings of the lowest degree modes 
(/ = 1 — 4). The data sets, available for this purpose are obtained by the ground- 
based networks BiSON IT, IRIS [61] and GONG ^45] and from the GOLF [79] 
instrument on SOHO. Unfortunately these sets of data are not in mutual agree- 
ment and give conflicting results of inversion in the core, as it is shown in Fig. 
[TBI taken from Di Mauro et al. [33j. Here, the radial spatial resolution, for clarity 
not drawn in the figure, is fixed at Ar = 0.1 Rq. The independent sets of obser- 
vations obtained by IRIS, GONG, GOLF lead to the conclusion that the Sun's 
core is in a state of rotation slightly faster than that observed at the surface in 
contradiction with the BiSON's data inversion which indicates a central angular 
velocity even slower than the surface polar angular velocity, as it was recently 
confirmed by Chaplin et al. [13] . Thus, the kinematics in the core remains largely 
uncertain, with a disagreement that might derive from the different data analysis 
procedures employed. 



30 Maria Pia Di Mauro 




Fig. 16. Rotation of the Sun's core as deduced by inversion of tlie BISON (filled 
circles), IRIS (filled triangles), GONG (starred symbols) and GOLF (filled squares) 
sets of lowest degree splittings (/ = 1 — 4), all combined with MDI higher degree data 
set. The radial spatial resolution of each radial point is fixed at Ar = 0.1 Rq 



7.4 Helioseismic determination of the solar angular momentum and 
quadrupole moment 

The present angular momentum of the Sun 3, can be deduced from the internal 
rotational behaviour derived from helioseismological data, by integrating the 
following ,33; : 

9=/ r^dMr (1- n)f2{r,n)dn = - f2i{rydMr , (68) 

Jo Jo •J Ja 

where f2i{r) is determined by helioseismic inversion of the Oi-splitting coefficient, 
from Eq. If we assume the angular velocity behaviour shown in Fig.[T31 the 
integration of Eq. §Sil leads to 9 = (1.96 ± 0.05) • lO*'^ g cm^ sec'^ [53]. This 
value is in agreement within errors with the one obtained by Pijpers in |72| . 

Another quantity of particular interest is the gravitational quadrupole mo- 
ment J2 of the Sun, which can be deduced, according to Pijpers [72| by evaluating 
the two-dimensional integral: 

J2 = ^ °drj J^{r,^i)n^{r,^l)df^, (69) 

where J-{r, fj,) is the two-dimensional kernel which depends on the physical quan- 
tities of the reference model and on some more general assumptions on the 
physics of the Sun. 
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The value of J2 obtained by Pijpers [TS] is J2 = (2.23 ± 0.09) x lO^'^. This 
result is totally consistent with the one obtained by Paterno et al. [69] with a 
different approach based on both the measurement of solar oblateness and the 
angular velocity profile deduced by inversion of splittings. 

8 Seismology of the fine structure: solar asphericities 

The asymmetric part of the fine structure in the p-mode spectrum (Eq. I63p of 
solar oscillations varies in a systematic way through the solar cycle [58] , [59] , [65] . 

It is evident that the changes are associated with the surface temperature 
bands reported by Kuhn et al. [60]. Also, Woodard & Libbrecht [95] found a 
strong correlation between oscillation frequency changes and solar surface mag- 
netic variations from monthly averages of their data. The origin of this behaviour, 
as well as the temporal variation of the frequencies is still ambiguous, but it ap- 
pears clear that all these changes are consistent with a near surface perturbation. 

The even-order splitting coefficients a2k.,i,n, seen in Eq. (|63l) can be fitted to 
the following formula obtained by Dziembowski & Goode [35] : 




where a2k,i,n;Tot represents the effect of centrifugal distortion which can be cal- 
culated following the treatment of Dziembowski & Goode S6.] /;^„ is a measure 
of the modal inertia; Ck,i is a constant which depends on the degree and on k, 
and 7fc is the asphericity coefficients which is directly related to the distortion 
described by the P2k{fJ-) Legendre polymonial. The P2{fJ-) term corresponds to a 
quadrupolar distortion (the oblateness), while PiifJ.) is the hexadecapole shape 
term and so on. 

Dziembowski et al. in [37] have studied the behaviour of the even splitting 
coefficients a2k,i.n of p modes, obtained by observations covering almost all the 
period during the past 11 years cycle. In Fig. [iTl taken from [37], there are 
shown the mean values of the jk coefficients as obtained by observations from 
various instruments, which include BBSO for the period 1986-1990, LOWL for 
the year 1994 and SOHO/MDI for the period 1996-1997. The variation of the 
asphericity coefficients is compared in Fig. [T7] with the monthly averages of 
smoothed sunspot numbers. Clearly, the BBSO data of 1988 and 1989 give the 
largest magnitudes of 71, 72 and 73, and this corresponds to the first half of the 
previous sunspot maximum. In years of high activity all three coefficients are 
substantial and change rapidly, while during low magnetic activity their value is 
roughly zero. In particular the asphericity appears more pronounced in period of 
high activity when it happens that P2ifJ') and P^ifi) distortions decrease while 
the PeifJ') distortion increases. This can be translated in the fact that the Sun 
assumes a shape which varies from simply oblate to complicated asphericity 
according to the magnetic cycle. 

This interesting conclusion has been confirmed by Howe et al. which 
analyzed data obtained by the GONG network during the period 1995-1998. 




a.2k,l,n — 



0-2k.l,n;iot + Ck,l— , 



(70) 
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Fig. 17. The lower panel is a smoothed monthly average of the sunspot number 
covering the time since the 1986 activity minimum. The upper panel is a combination 
of 7's derived from BBSO (1986-90), LOWL (1994) and SOHO/MDI (1996-97) data. 
The errors in the 7's are smaller than the symbols used to represent their values [37) 



They also observed that the temporal variation of the a2fe-coefhcients is strongly 
correlated with the latitudinal distribution of the surface magnetic activity. 

The behaviour of the 7's as function of the frequency [37] yields, also, informa- 
tion about the sources of the solar distortion reflected in the even-a coefficients. 
It is well-known that p modes sample the region just above their inner turning 
points. Recently, Dziembowski et al. 'SiS'i found a significant aspherical distortion 
in the layer located at a depth ranging between 25 and 100 Mm. The perturba- 
tions seems to arise from a relative temperature increase of about 1.2 x 10^^ or 
from a magnetic perturbation, with (B^) ~ (60 KG)^. 

9 Concluding remarks 

Helioseismology, through the very accurate identification of oscillation frequen- 
cies of acoustic and fundamental modes, has clearly demonstrated that the stan- 
dard solar models reproduce the behaviour of the Sun with remarkably accuracy, 
consistent within 1 %. 

Despite such overall success, this discipline has not yet exhausted its re- 
sources, since helioseismic results clearly suggest further refinements of the solar 
models. 

The detailed structure of the convective zone and of the near-surface region 
is quite uncertain, since there remains substantial ambiguity associated with 
modelling the convective fiux, taking into account the non-adiabatic effects, ex- 
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plaining the excitation and damping of the solar oscillations and defining an 
appropriate equation of state to describe the thermodynamic properties of the 
solar structure. 

The attempts to restore the solar core conditions, up to now, have been 
contradictory too. In fact p modes (as opposed to gravity modes, g modes) are 
not very sensitive to the core of the Sun. This indicates the necessity of using 
more accurate low degree p-mode data and to continue to investigate for the 
presence of g modes. 

In addition, there is still much work ahead in getting a detailed understanding 
of the Sun's rotation. Some rotational features like, for example, the temporal 
changes which occur near the base of the convective envelope have not been yet 
explained. 

Finally, by studying the connection between the seismic and the global char- 
acteristics of the Sun, the challenge is to find the reason for the correlation 
between the variation of the Sun's shape and the magnetic solar cycle. 

Ever more precise helioseismic observations from ground and space can help 
us to reconstruct the complete picture of the Sun and, finally, to solve the most 
discussed open questions in solar physics such as the solar neutrino problem, the 
history of the Sun's angular momentum, and the solar cycle generation mecha- 
nism, through the interaction of the convective motions with the rotation inside 
the Sun. 

Recently, a new window has been opened on the astrophysics research: the 
possibility to study and to understand the behaviour of other stars by applying 
the tools and the techniques well developed and used in helioseismology. In fact, 
the success of helioseismology has spurred investigators to extend this diagnostic 
to other stars which may show multi-mode pulsations. Up to now, the seismo- 
logical study of pulsating stars, known as Asteroseismology, has been hindered 
by the problem of mode identification since the oscillation amplitudes observed 
on the Sun (a few parts per million in flux) are too small to be detected in other 
stars with ground-based telescopes. To reach the required sensitivity and fre- 
quency resolution, several space experiments, MONS [54], COROT [S], MOST 
[66] . will soon be devoted to the measurements of stellar oscillations. Thus, it is 
evident that asteroseismology represents the successive step in the evolution of 
the helioseismology research. 
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